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Abstract. Isospin-breaking effects, in particular those associated with electroweak- 
I , penguin contributions and ir°-r], rf mixing, have long been known to affect the Standard 

' Model prediction of e'/e in a significant manner. We have found an heretofore uncon- 

sidered isospin-violating effect of importance; namely, the u-d quark mass difference 
can spawn |A/| = 3/2 components in the matrix elements of the gluonic penguin op- 
erators. Using chiral perturbation theory and the factorization approximation for the 
hadronic matrix elements, we find within a specific model for the low-energy constants 
that we can readily accommodate an increase in e'/e by a factor of two. 

o ; 
^ : 

O*'. INTRODUCTION 

The recent measurement of a non-zero value of Re (e'/e) [1] establishes the ex- 
istence of CP violation in direct decay and thus provides an important first check 
| of the mechanism of CP violation in the Standard Model (SM). Nevertheless, the 

world average which emerges is Re (e'/e) = (19.3±2.4)-10 -4 [2], which is larger than 
the "central" SM prediction of 7.0- 10~ 4 [3,4] by nearly a factor of three. This com- 
pels us to scrutinize the SM prediction in further detail: we study isospin-violating 
effects arising from the u-d quark mass difference. 

Isospin violation plays an important role in the analysis of e'/e, for the latter is 
predicated by the difference of the imaginary to real part ratios in the |AJ| = 1/2 
and | A/ 1 = 3/2 K —>■ tttt amplitudes. The differing charges of the u and d quarks 
engender |A/| =3/2 electroweak penguin contributions, whereas 7r -i],ri' mixing, 
driven by the u-d quark mass difference, modifies the relative contribution of the 
| AI | = 1/2 and |A/| =3/2 amplitudes in a significant way. 

Here we describe isospin-breaking effects in the matrix elements of the gluonic 
penguin operators [5], such as Qq. These operators have always been thought to 
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induce exclusively |AJ| = 1/2 transitions, but this is true only in the limit of isospin 
symmetry. The difference in the up and down quark masses effectively distinguishes 
the interaction of gluons with up and down quarks, so that the (tt7t\Qq\K) matrix 
element possesses a |A7| = 3/2 component as well [6]. 

Let us begin by showing why the numerical prediction of e'/e is sensitive to the 
presence of isospin violation. The value of e'/e is inferred from a ratio of ratios, 
namely 
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In the isosp in-perfect limit, the two independent amplitudes present in K — > Tin 
decay are distinguished by the isospin of the final-state pions, namely Ai = A(K — > 
(7T7r)/) with / = 0, 2. e'/e can thus be written 



e 
e 



r e(i-n), 



with 



ReA 2 
ReA 



£ = 



Im A 
ReA 



n = 



Im A 2 
culm A 



(3) 



(4) 



In standard practice, ui ~ 1/22 and Re A are taken from experiment whereas ImAj 
is computed using the operator-product expansion [3,4], that is, via 
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H eS (\AS\ = l)=4^=v: s V ud J2C t ^)Q^) + h.c. 
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The numerical value of e'/e is driven by the matrix elements of the QCD pen- 
guin operator Q$ and the electroweak penguin operator Q s [7]. Writing (Qi)i as 



((tttO^IX) = B 



((J+l)/2) 



((tttt)/^^)^), where ' 'vac" indicates the use of the 



vacuum saturation approximation, one recovers the schematic formula [3] 
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Using B£ /2) = 1.0, B£''> = 0.8, and Q v+ri , = 0.25 yields the "central" SM value 
of e'/e ~ 7.0 • 10~ 4 [3], roughly a factor of three smaller than the measured value. 
Larger estimates of B^ 2 \ and hence of e'/e , exist [8-10]; we investigate sources 



of J\ +r) '. Note that under Q v+V > — > — Q v + V ', Re(e'/e) — > 2.2 Re(e'/e). Were 
l-B^^I ^ l-^8 3 ^ 2 ' ) |) flipping the sign of Q v+V ' would increase e'/e by a factor of 1.7. 
Let us consider possible sources of fl v +ri'- We replace Q v +-n' by Qm, where 



and lmA P (K° -> tttt) = (ImA P (ir° -> tt+tt-) + ImA P (K° -> 7r°7r°))/2. u A p 



denotes an amplitude induced by (8l, 1r) (e.g., Qq) operators - - the empirical 
|A/| = 1/2 rule suggests such operators dominate the isospin- violating effects. 
Qib vanishes in the absence of isospin violation, i.e., if m u = m d , e u = e&. It can 
be generated by both strong-interaction and electromagnetic effects, mediated by 
m d 7^ m u and e u ^ e d [11], respectively. We focus on m d ^ m u effects. The latter in- 
clude ir°-rj,rf mixing [12-14]; in 0(p 2 , 1/N C ) this yields Q v+V > = 0.25 ±0.05 [13,14], 
used in the analysis of Ref. [3]. However, m u ^ m d effects can also spawn a 
| AI | = 3/2 component in the matrix elements of the gluonic penguin operators [6], 
as illustrated in Fig. 1. We turn to a chiral Lagrangian analysis in order to estimate 
the size of this effect [5]. 



The weak chiral Lagrangian for K — > tttt decay is written in terms of the unitary 
matrix U = exp(i<f>/f) and the function x, both of which transform as U — > RUL^ 
under the chiral group SU(3)l x SU(3)r. The function represents the octet of 
pseudo-Goldstone bosons, i.e., = Z) a =i ... 8 ^a<Pa- In the absence of external fields, 
X = 2BqM with M = diag(m u , m s ) and Bq oc (qq). The leading-order, 0(p 2 ), 
weak chiral Lagrangian contains no mass-dependent terms [15], so that ^ m u 
effects in the hadronization of the gluonic penguin operators first appear in 0(p 4 ). 
This is illustrated in Fig. 2. 



FIGURE 1. Quark line diagram illustrating the "strong penguin" s — > dqq transition in the 
Standard Model. Note that qq € uu,dd; in the isospin-perfect limit, m u = and only 
|A/| — 1/2 transitions are generated. If m u ^ nid, the K —> tttt matrix element associated 
with this operator contains a |A7| = 3/2 component as well. 





CHIRAL LAGRANGIAN ANALYSIS 
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Let us enumerate the possible isospin- violating effects which occur in 0{m d —m u ) 
and 0(p A ) : 

i) 7r°-i] mixing realized from the 0(p 2 ) strong chiral Lagrangian, in concert with 
the 0(p 2 ) weak chiral Lagrangian, computed to one-loop order. 

ii) 7c -r] mixing, realized from the 0(p 2 ) strong chiral Lagrangian, combined with 
the isospin-conserving vertices of the 0(p 4 ) weak chiral Lagrangian. 

iii) TT°-r] mixing as realized from the strong chiral Lagrangian in 0(p A ), combined 
with the 0(p 2 ) weak chiral Lagrangian. The 7r -?/ mixing effects included in 
Refs. [13,14] ape this effect. 

iv) Isospin violation in the vertices of the 0(p 4 ) weak chiral Lagrangian. This 
serves as our focus here, for it contains the qualitatively new effects we argue. 

We use the octet terms in the 0(p A ), CP-odd weak chiral Lagrangian of Ref. [16]. 
Collecting the %-dependent terms as per iv), working to 0{m d — m u ) ) and dropping 
terms suppressed by M 2 jM\, we find 

_ 2^2 M 2 K0 B (m d -m u )~ 0.12GeV 2 ~ 

iZp = k hi ~ hi o 

3lu M 2 k0 - M 2 c 2 - c 2 1 ; 

with E~ = 2E1-2E3 -4£?4 -E w -E u -AE^-E^. Note that c 2 is the low-energy 
constant associated with the 0(p 2 ), (8l,1r) weak chiral Lagrangian [16]. As per 
ii), 7r°-?7 mixing in 0(p 2 ) also enters when combined with the isospin-conserving 
vertices of the 0(p 4 ) weak chiral Lagrangian. Including the x-dependent octet 
terms, we find 
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FIGURE 2. Isospin violation in K — > rnr decays. The square box represents the weak \AS\ — 1 
transition at low energies, whereas the "x" represents the presence of 7^ m u effects. Mass 
effects do not occur in the weak transition in leading order in chiral perturbation theory, so that 
only the left-hand diagram occurs in 0{p 2 ) and in 0(rrid — m u ) — 7r°-?7 mixing, mediated by 
md — m u effects in the strong chiral Lagrangian, can occur. In 0(p 4 ) both diagrams are present; 
the new effect we discuss is associated with the right-hand diagram. 



with E- +7)l = -2{Es + El - £ 5 ~) + (E{ - E^)/ 2 - 2E U + £f 4 + 3£f 5 /2 — so that 
no manifest cancellation with the terms of Eq. (8) occurs. 

The low-energy constants E^ are unknown, so that we turn to the factorization 
approximation to proceed. The construction relevant to (8^, 1r) transitions in 
K° — ► 7T7T decay is [17] 

£ p= ~^l V :sVudC 6 (-8(s L q R )(q R d L )^ + h.c 
- %V: s V ud C^2Bl^^ + h.c. , 

where £ str is the strong chiral Lagrangian. To generate terms of 0{p i ) in Cp 
requires terms of both 0{p A ) [18] and 0(p 6 ) [19] in £ str . Unfortunately, the low- 
energy constants of the latter are also unknown; the use of "resonance saturation" 
allows us to estimate some of them. We explicitly consider the scalar nonet of res- 
onances as per Ref . [20] . An example of the manner in which the scalar resonances 
can generate contributions to the E~ is illustrated in Fig. 3. Integrating out the 
scalar resonances for p 2 <C M|, we find two terms which contribute to the scalar 
densities in the bosonization of Qq [20], 

A 6 ) — ^ mC m /^3\ | c dC m d m 2 2 

Ls ~ 2Mj {X+} + ~MT {X+ } ' 1 } 

yielding contributions to E{ and Ei in terms of d m , c m , Cd, and Ms- The pa- 
rameter d m is ill-known; we find d m ~ —2.4 (—0.76). The sign of d m and thus 
of fip in our model results from the mass of the lowest-lying strange scalar being 
greater than that of the lowest-lying isovector scalar. As per our earlier classifica- 
tion, ftjg = Qjg'* + ftjjlp + Q$' iU + Q<$' m , so that with d m = -2.4 (-0.76), we 
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FIGURE 3. A contribution to the right-hand diagram of Fig. 2 in 0(rrid — m v ), estimated 
in the factorization approximation with explicit, scalar-resonance degrees of freedom. The "(g)" 
represents a bosonized current; the open parentheses denote contributions from the vertices of 
the 0(p 6 ) model La grangian of Ref. [20] . In this case the isospin- violating contribution is driven 
by a -fo mixing. 



have 4 " = —0.79 (—0.21). Estimating f% using the x-dependent E i yields 
^ib'" = —0.12 (—0.03). fijg'™ has been partially determined through the inclusion 
of 7T°-V mixing in Q v+T] , = 0.25 ± 0.05 [13,14]. Using the result fig + fig'^ = 
0.16 ± 0.03 [21] and neglecting fig' 1 , as the ill-known E~ do not warrant such a 
calculation, we estimate, finally, that fins = fins + fins ~ —0.05 — > —0.78. For 
reference, note that fig ~ 0.13. The large value of fig is driven by the numerical 
prefactor of Eqs. (8,9) — the contributions in fig* are "naturally" of the same size as 

(2) 

JIjb ; . Thus we find a very large correction to the value of £l v+v i = 0.25±0.05, used in 
"central value" of e'/e . The large negative change in fi IB found in 0{p A ) generates 
a substantial increase in e'/e . 

The fi ffi we calculate impacts e'/e in a significant manner. Our estimate of fii B 
from the specific m& ^ m u effects we consider ranges from —0.05 — > —0.78; this 
range exceeds the central value, £l v +r/ — 0.25 ± 0.05, used in earlier analyses and 
reflects a variation in e'/e of more than a factor of two. The presence of unknown 
low-energy constants implies that we lack a reliable way to calculate the effects we 
consider. Such limitations, however, underscore the need for a larger uncertainty 
in the Standard Model prediction of e'/e . 
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